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1  Introduction 


Consider  k  independent  logistic  populations  Hi, . . . ,  Ilfc  with  unknown  means  (9i, . . . ,  4- 
Let  0[i]  <  . . .  <  0[k]  denote  the  ordered  values  of  the  parameters  It  is  assumed 

that  the  exact  pairing  between  the  ordered  and  the  unordered  parameters  is  unknown. 
A  population  Ilj  with  9i  =  9[k]  is  called  the  best  among  the  k  underlying  populations. 
In  many  pratical  situations,  we  may  not  only  be  interested  in  the  selection  of  the  best 
population,  but  also  require  the  selected  population  to  be  good  enough  compared  with 
a  given  control.  The  problem  of  selecting  the  best  population  has  been  studied  by  many 
researchers.  Gupta  and  Panchapakesan  (1996)  provided  a  comprehensive  review  of  the 
development  in  this  area.  It  should  be  pointed  out  that  the  logistic  distribution  serves  as 
a  statistical  model  in  many  practical  situations,  see,  for  example,  Balakrishnan  (1992). 
The  statistical  selection  problem  for  logistic  populations  has  been  studied  in  Gupta  and 
Han  (1991, 1992),  among  others. 

In  this  paper,  we  employ  the  empirical  Bayes  approach  to  select  the  best  logistic 
population  provided  it  is  also  as  good  as  a  given  control.  We  describe  the  formulation  of 
the  selection  problem  and  derive  a  Bayes  selection  procedure  in  Section  2.  In  Section  3, 
we  construct  an  empirical  Bayes  selection  procedure.  Then  we  investigate  the  asymptotic 
optimality  of  the  proposed  empirical  Bayes  selection  procedure  in  Section  4.  A  simulation 
study  is  carried  out  to  investigate  the  performance  of  the  proposed  selection  procedure 
in  Section  5. 


2  Formulation  of  the  Selection  Problem 

Let  Hi, . . . ,  life  be  k  independent  logistic  populations  with  unknown  means  9i,  ...,9k- 
Let  9[x]  <  ...  <  9[k]  denote  the  ordered  values  of  the  parameters  9x,. ..  ,9k.  It  is  assumed 
that  the  exact  pairing  between  the  ordered  and  the  unordered  parameters  is  unknown. 
A  population  TTj  with  9i  —  9[k]  is  considered  as  the  best  population.  For  a  given  control 
9q,  population  tt^  is  defined  to  be  good  if  the  corresponding  9i>  9q,  and  bad  otherwise. 
Our  goal  is  to  select  a  population  which  is  the  best  among  the  k  populations  and  also 
good  compared  with  the  standard  ^o-  If  there  is  no  such  treatment,  we  select  none. 

Let  =  (^1, . . . ,  0fe)}  be  the  parameter  space.  Let  a  =  (uq,  ■  ■  ■ ,  Qfc)  he  an  action, 

k 

where  a,  =  0, 1; i  =  0, 1, . . . , /c  and  Y,  Oi  =  1-  For  each  i  =  1, . . .  ,k,ai  =  1  means  that 

1=0 

population  TTj  is  selected  as  the  best  and  also  considered  to  be  good  compared  with  9q. 
Co  =  1  means  that  all  the  k  populations  are  excluded  as  bad  and  none  is  selected.  We 
consider  the  loss  function 


1 


k 


L{9,  a)  =  max(%],  Bq)  -  ^  aiOi. 

i=0 


It  is  the  absolute  error  loss. 

For  each  i  =  1,...  ,fc,  let  be  a  sample  of  size  M  from  the  logistic 

population  Ilj  =  which  has  unknown  mean  9i  and  unknown  variance  (7r^cr?)/3, 

that  is,  the  conditional  density  distribution  of  Xij  given  9i  and  cr?  is 


1 

CTj  (1  +  ’ 


— oo  <  Xi  <  oo. 


(1) 


Since  logistic  distribution  is  symmetric  about  its  mean,  the  mean  and  the  median  of 
a  logisitc  population  distribution  are  identical.  For  convenience,  suppose  M  is  an  odd 
number,  and  we  denote  M  =  2s  +  1.  We  also  assume  that  the  unknown  population 
median  (and  also  the  mean)  9i  has  a  normal  N{fXi,Tf)  prior  distribution  with  unknown 
parameters  {ni,  r?).  The  random  variables  9i,  ...,9  k  are  assumed  to  be  mutually  inde¬ 
pendent.  Define  Xi  to  be  the  median  of  {Xu, . . .,  XiM},  i  =  1,. . .  ,k.  Let  fi{xi\9i,  of) 
and  hi{9i\iJii,  Ti)  be  the  conditional  distributions  of  Xi  given  {9i,  of)  and  9i  given  (^j,  rf), 
respectively.  We  have,  for  i  =  1, . . . ,  A:, 


ji{Xi\9i,(jl)  = 


(2s  +  l)!l 

(s!)^  CTj  (1  -l-  e~(®*~^<)M)2«+2  ’ 


—  oo  <  Xi  <  oo. 


(2) 


From  (2)  we  see  that  the  density  function  fi{xi\9i,af)  is  symmetric  about  9i  given  9i, 
therefore. 


EXi  =  E{E{Xi\9i))  =  E9i  =  iXi. 
The  posterior  density  of  9i  given  Xi  =  Xi  is  proportional  to 


(1  -1- e-(a:<-5i)/o-j^2s+2  ^ 


— OO  <  <  OO. 


(4) 


Let  X=(Xi, . .  .,Xk)  and  X  be  the  sample  space  generated  by  X.  A  selection  pro¬ 
cedure  d=(do,  ■  ■■,dk)  is  a  mapping  defined  on  the  sample  space  X.  For  every  x  E  X, 
di{x),i  =  1, ...  ,k,  is  the  probability  of  selecting  population  Ilj  as  the  best  among  the  k 
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populations  and  also  good  compared  with  the  given  control  0o,  do(x)  is  the  probability 
of  excluding  all  k  populations  as  bad  and  selecting  none.  Also,  J2i=odi{z)  =  !>  for  all 
xGAf. 

Under  the  absolute  error  loss,  the  posterior  median  is  the  Bayes  estimator  of  9i.  We 
denote  (pi{xi)  to  be  the  posterior  median  of  9i  given  Xi  =  Xi,  i  =  I, . . .  ,k. 

Under  the  preceding  statistical  model,  the  Bayes  risk  of  the  selection  procedure  d  is 
denoted  by  R{d).  We  have 


R{d)  =  -  /  [Y,di{x)(pi(xi)]f{x)d{x)  +  C, 

i=0 


where 

C  =  /n  max(0[fc],  9o)dH{9), 

H{9):  the  joint  distribution  of  9  =  {9i, . . . ,  9k), 

=  !Rfiixi\9i,a^)hi(9i\ni,Tl)d9i, 
f{x)=IllJi{Xi), 

<Po(xq)  =  9o. 

For  each  x  G  A,  let  /(x)  =  {i\ipi{xi)  =  mRX(pj{xj),i  ~  0,1,..., A:},  and  i*  = 

min{i|i  €  I(x)}.  Then  a  Bayes  selection  procedure  d^{x)  =  (d^(x), . .  .,dk(x))  is  given 
as  follows: 


df.(x)  =  1, 
df(x}  =  0,  for 


3  The  empirical  Bayes  Framework 

The  Bayes  selection  procedure  d^(x)  defined  in  Section  2  depends  on  the  unknown 
parameters  i  =  l,...,k  and  the  specific  form  of  ^i{xi).  Since  the  parameters 

and  the  specific  form  of  (pi{xi)  are  both  unknown,  it  is  impossible  to  implement  the 
Bayes  selection  procedure  for  the  selection  problem  in  practice.  In  the  empirical  Bayes 
framework,  it  is  generally  assumed  that  there  are  some  past  observations  when  the 
present  selection  is  to  be  made.  At  time  I  =  1, . . . ,  n,  let  Xiji  be  the  j-th.  observation 
from  Hi,  that  is,  for  each  z  =  1, . . . ,  A:,  let 


3 


and 


(7) 


Xiji  ~  L{diu  O’- ),  j  = 


For  I  =  1,. .  .,n,  denote  Xi^i  to  be  the  median  of  {Xm, Xmi),  and 


(8) 


Xi{n)^-j^Xi^u 

1=1 


(9) 


SKn) 


1 

n  —  1 


^(X,J  -  Xi(n))^ 


l=l 


Then, 


(10) 


E{Xi,i)  =  E{E{Xi,i\9u))  =  E{9u)  =  Mi, 


(11) 


and 


Var{Xi,i)  =  Var{E{Xi,i\9u))+E{Var{Xi,i\9u)) 

=  Var{9u)  +  E{Var{Xi,i\9u)) 

=  Tf  +  E{Var{XM) 

<  oo.  (12) 

Denote  uf  =  Var(Xi^i).  Since  (Xa,...,Xin)  are  i.i.d.,  by  the  strong  law  of  large 
numbers,  we  know  that  as  n  — )■  oo, 

f  -^i(’^)  ^  Mi, 

I  (n)  — u?,  a.s. 

To  derive  the  empirical  Bayes  selection  procedure,  we  first  consider  the  following 
lemmas.  The  following  lemma  is  from  Serfling  (1980). 
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Lemma  3.1  Let  {Fi,  1  <  i  <  mj-  be  m  i.i.d.  random  observations  from  continuous 
distribution  function  F;  also  let  |  and  ^  be  the  medians  of  {Yi,l  <  i  <  m}  and  F, 
respectively.  Then,  for  any  e  >  0, 


- -^1  >  4  < 


(14) 


where  5^  =  min{F(^  +  e)  “  2’  2  “  “  ^)}- 

Put  a'  =  mini<j<fc  (T,,  a*  =  maxi<i<fc  cTj.  Xu, XiM  are  i.i.d.  from  L{6i,  cr?),  which 
has  the  following  cumulative  distribution  function 


m  = 


1 

1  -j-  Q~{U—9i)l<^i 


—  OO  <  ti  <  CO, 


(15) 


and  for  0  <  e  <  a', 


jC/o-i  _  1 


2(e«M  +  1)  “  2(e  +  l)cr*‘ 


> 


(16) 


Given  9i,  6i  and  Xi  are  the  population  median  and  sample  median  respectively,  we 
have,  from  Lemma  3.1, 


-(23  +  l)t^ 

P{\Xi  -  9i\  >  e}  <  (17) 


For  any  0  <  e  <  a',  denote  =  {x  6  :  \xi  —  9i\  <  e}.  We  show  that  the 

conditional  density  of  Xi  given  9i  and  erf  is  approximately  N{9i,  yfycrf)  as  s  — f  oo. 

Prom  (2),  the  conditional  density  of  Xi  given  9i  and  erf  is 


fiixi\9i,ai)  = 


(2s  +  l)!l  (e-(®i-»dM)s+i 

(s!)^  CTj  (1  +  e“(^*“®d/<ri)2s+2 

(2s  +  l)!  1 _ 1 _ 

(s!)2  (Ti  (2  +  4-  eixi-9i)/aiy+i ' 


By  Stirling’s  formula,  when  s  is  large  enough. 


(18) 
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(19) 


(2s  + 1)! 
(s!)2 


2(2^+§)  - - - 

Ri  — -^=^Vs  +  1. 


Also  choosing  e  =  4-  0  to  be  a  sequence  of  fixed  numbers  which  tend  to  0  as  s  — )•  oo, 

by  Taylor’s  polynomial  expansion,  we  have 


log(2  +  log4  +  (20) 

4  (Tj 

on  Si-  When  s  ->  oo,  from  (17), 

-(23+1)£2 

P{x  ^  5i}  <  — )■  0.  (21) 

Therefore,  we  see  that  as  s  — >■  oo, 


fi{xi\6i,a'l) 


.iillfizM; 
4  \r2 


^/^^J2/s  +  l<^i 


(22) 


that  is,  /i(a;i|^i,  cr?)  is  approximately  N(di,-^crf). 

From  above,  we  can  see  that  for  sufficiently  large  s,  the  conditional  density  of  Xi^i 
is  approximately  N(di,  j^crf),  given  Oj  and  cTj-.  Since  the  prior  distribution  of  0i  is 
N{Hi,Ti),  the  unconditional  density  of  Xi^i  is  approximately  N{ni^Tf  +  -^crf). 

For  each  population  flj,  let  Wi{n)  be  the  measure  of  the  overall  sample  variation  for 
the  past  observations.  That  is. 


I  =  (wiiF  Eiii  a.Wii  - 

Then  we  define,  for  i  =  1, . . . ,  A:, 

Ai  ~  A^((7l), 

^ 

vf  =  5?(n), 

.  f?  =  max(i>?  -  0). 


(23) 


(24) 
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and 


{Xirf  +  ^fti)li>h  if  - 


l^ii 

I  0q{xq)  =  do. 


if  A 


.•',2 


s+1"* 


a}  <  0, 


(25) 


Then  for  each  x  €  X,  let  /(x)  =  {i\(pi{xi)  =  m!^^(pj{xj),i  =  0,1,..., A:},  and 

i*  =  min{?|i  €  /(x)}-  We  propose  the  following  empirical  Bayes  selection  procedure 
^(n,s)(^)  _  (x), . . . ,  as  follows: 


I  =  1, 

jdM  =  0,  iov  j^i*. 


4  Performance  of  the  proposed  selection  procedure 

Consider  the  empirical  Bayes  selection  procedure  di”’*^(x)  constructed  in  Section  3. 
Let  be  the  conditional  Bayes  risk  given  the  past  observations  {Xiji,i  — 

1,. . .  ,k;j  =  1, . . . ,  M;  and  I  =  1, . . . ,  n}  and  ER{dX^'^\x))  the  Bayes  risk  of  the  em¬ 
pirical  Bayes  selection  procedure  respectively,  where  E  is  the  expectation  taken  with 
respect  to  the  past  observations  {Xiji}.  From  (5), 


i2(d(”’"^(x))  =  -  /  E  df'^\x)>pi{xi)]f{x)d[x)  +  C.  (27) 

i=0 

Note  that  i?(d("'’^)(x))  —  R(d^(x))  >  0,  since  d^(x)  is  the  Bayes  selection  procedure. 
Therefore,  E(R(d^’^’^^(x))  —  R(d^(x)))  >  0.  We  use  the  nonnegative  difference  regret 
risk  E{R{d^^’^'>{x))  -  i?(d^(x)))  >  0  as  a  measure  of  the  performance  of  the  selection 
procedure  d^'^’^^x). 

We  first  state  some  facts  about  (pi{xi),  the  posterior  median  of  dj  given  Xi  =  Xi  and 
From  the  definition  of  ipi{xi),  we  can  see  that  (pi{xi)  is  between  Xi  and  fXi.  Besides, 

Lemma  4.1  When  s  is  large  enough,  for  1  <  i  <  A:, 


\fPi{xi)  -Xi\<  2(Xi 


(28) 
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Proof.  We  only  prove  (pi{xi)  <Xi  +  here.  The  proof  of  ^i{xi)  >Xi  — 

is  similar.  To  prove  ^i{xi)  <  Xj  +  2<rjy  it  suffices  to  show  that 


roo 

/  fi{xi\ei,al)hi{ei\nu  rf)dei 

Jxi+2cTi\/ 

/•°°  (2s  +  l)!l  1 

Jxi+2<Ti-y^^  (s!)2  (Tj  V^Ti  (1  + 

1  (2s +  1)!/  e® 

(5!)2  l(l  +  e'^)2j 


de-^0, 


(29) 


as  s  — )■  oo.  We  first  show 


t(e,s)  :=  ^  — >0  ass^oo,  (30) 

uniformly  for  6  >  2\f^^.  Obviously  it  is  enough  to  consider  the  case  oi  9  — 

since  t{9,  s)  is  decreasing  on  0  >  0.  When  9  =  2-^^^  and  s  is  large  enough,  by  Taylor’s 
formula, 


log(l  +  e®)  —  log  2-\-  -9  +  -9^  +  o{9^) ,  (31) 

and  by  (19),  when  s  is  large  enough, 

<  2(s  +  l)log2  +  ^log(s  +  l).  (32) 

From  (31)  and  (32),  we  obtain  that 


logt(0,s) 

=  (s  +  1)  [^  -  2  log(l  +  e®)]  +  log 

S  4-  1  1 

<  -2(s  +  1)  log  2 - —9“^  +  2(s  +  1)  log  2  +  -  log  s  +  o{s9‘^) 

=  -(^-^  -  i)logs  +  o(logs) — >--oo. 


(33) 
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as  5  oo.  Therefore,  (30)  is  proved,  from  which  we  can  immediately  see  that  (29)  holds 
true.  It  completes  the  proof  of  Lemma  4.1. 

The  next  lemma  is  well  known  and  can  be  found  in  Baum  and  Katz  (1965). 

Lemma  4.2  Let  Xi, . . . ,  be  i.i.d.  random  variables  with  mean  0.  Suppose  for  a  >  1, 
E\Xi\°‘  <  oo,  for  i  =  1, . . . ,  n,  then  for  any  e  >  0, 


(34) 

i=l 

As  a  consequence  of  Lemma  4.2,  we  have 

Lemma  4.3  Let  Xi, . . .  ,Xnhe  independent  random  variables,  with  mean  EXi  =  //  and 
variance  VarATj  =  for  i  =  1, . . . ,  n.  Also  let  X  =  ^YlXi  and  ^  T,(Xi  —  X)^. 
Suppose  for  z  =  1, . . . ,  n  and  a  fixed  number  a  >  2,  .SjA’il"  <  oo,  then  for  any  e  >  0, 


=  (35) 

The  proof  of  Lemma  4.3  can  be  found  in  Gupta  and  Lin  (1997). 

Since  EXf  i  <  oo,  for  any  e  >  0,  by  Lemma  4.2, 


P{|A  -  «l  >  £}  =  *). 


(36) 


also  by  Lemma  4.3, 


PM  -  ^f\  >  4  -  o{n-^). 


Similarly,  we  have  for  any  e  >  0, 


(37) 


P{\S-I  -  <^11  >  £}  =  »("-')■ 


(38) 


When  s  is  large  enough,  uf  —  -^crf  >  0.  Therefore,  from  (37)  and  (38),  when  s  is 
sufficiently  large. 


s  +  l"^* 


<  0}  =  o{n~^). 


(39) 
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Besides,  rf  =  uf  —  E{Var{Xi^i\9i))  by  (12)  and 


/OO 

(xii  -  ei) 

■00 


We  have 

Lemma  4.4 


2  (25  +  1)!  1 

(s!)^  cTj  (1  +  e~(®*'“^dM)25+2 

•/-»  (s!)2  V(l  +  e-)2' 


(40) 


Proof. 


/_ 


=  2f 

Jo 


(51)2  V 

(1  +  e^)2i 

2(25  +  1)1 

/  e"  \ 

'  (^0^ 

V(l  +  e’')2i 

,(2s  +  l)!/  e*  ' 

s  J 


s+l 

dx 

s+l 

dx 


rV^-^  ^  ^  roo  ^^{2s  +  l)\ 

Vo  is 


(5!)2  V(l  +  e*)2 


s+l 


dx 


:=  Ti+Ta+Ts. 

By  Stirling’s  formula,  when  s  is  large  enough, 


(41) 


(42) 


Ti  =  2 


(25  +  1)!  6“ 


/' 


X 


(1  +  e®)2 


a+l 


dx 


<  2  ■  22(^+')v7+I  •  2-2(^+i)  ["^x^dx 

Jo 


< 


=  01 


gl0gS^3/2 


flog  5 


Using  the  same  approach  as  in  the  proof  of  Lemma  4.1,  we  have 

,(25  +  1)1  y  2 


(s!)2  iyil^  V(l  +  eVV 


s+l 


(l  +  e“)2. 


dx 


(43) 
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<  2 


=  0 


(2s +  1)!/  e 


\«+i  /-s 

/  , _ X  dx 


(s!)2  V(i  +  e^iT)2/  JV^- 


flogs 


(44) 


Moreover, 


=  2 


(2s +  1)!  /•“ 


r^G 


s+l 


dx 


(s!)2  Jz  V(l  +  e®)2 

(s!)^  Jz 
-  “(#)■ 


This  completes  the  proof  of  Lemma  4.4. 

Prom  Lemma  4.4,  we  observe  that  when  s  is  sufficiently  large. 


(45) 


E{Var{Xi,i\ei))  =  o 


(logs 


(46) 


and  therefore,  by  (37),  (39)  and  the  definition  of  ff ,  for  e  >  cy^^,  where  c  >  0, 


P{\n-r^\>^}  =  o{n-^) 


(47) 


and  furthermore. 


P{c>'tM  <  =  »(»-')■  (48) 

Next  we  investigate  the  rate  of  convergence  of  E{R{d^'^'^'>{^)—R{d^{x))).  Let  P„,s  be 
the  probability  measure  generated  by  the  past  observations  Xiju  i  =  1, . . .  A;;  j  =  1, . . .  M 
and  Z  =  1, . . .  ,n. 


P(P(d("’"nx))  -  i?(<^^(x))) 


hi*  =  3}{^i{xi)  -  (Pj{xj))f{x)dx 
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=  E  /  =  i,i*  =  0}(¥>i(a:i)  -  9o)fix)dx 

+  E  /  PnA^*  =  0,  f  =  iX^o  -  ^Axj))f{^)d-i 

j=\ 

+  E  E  /  PnA'^  =i,l*  =  j}{<PiiXi)  -  ipj{x^))f{x)dx 

i=l  j=\ 

k  ^ 

^  2^^  J ^Pn^s{\^i{^i)  ~  ¥^i(^i)|  ^  ~  ^o|}|*^i(^i)  9Q\fi[Xi)dXi 

+2  E  E  /  -  wfe)i  >  -  >pi(^i)\ 


:=  /1  +  /2. 


X  fi{xi)fj{xj)dxidxj 


(49) 


For  any  e  >  0,  and  i,j  —  l,...,k,  let 


—  {^i  ■  ^o|  ^  ^}) 

Xij  =  {{xi,Xj)  :  \^Pi{xi)  -  (Pj{xj)\  <  e}. 


(50) 


Then  we  have 


-^1  =  ^E  /  PnA\Vi{^i)  -  ^i{^i)\  >  -  ^o|}iV?i(a;i)  -  %\fi{Xi)dXi 

i-X 

^  /• 

+2E  /  ■Pn,XI<i^i(3^i)  -  <^i(^^)l  >  bi(^i)  “  ^0 1 } 1 95i (^^i )  “  ^o|/i(a^i)rfa;i 

i=l 

<  2E  /  efi{xi)dxi 

i=l 

k  n 

+2E  j^PnAWii'^i)  -  '^i(^i)l  >  ^}\Vi{^i)  -  9Q\fi{Xi)dXi.  (51) 


By  Lemma  4.1,  when  s  is  large  enough,  \^Pi{xi)  —  a;,]  <  From  now  on,  we 

always  set  e  =  Therefore,  for  sufficiently  large  s, 


Xi-Ool  <  \ipi{xi)  -  Xi\  +  \ipi{xi)  -9o\<2e 


on  Xi  and 
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JXi  J{\xi-0o\<2e} 


<  f 

J(\X4-' 


{|xi-«o|<26}  a/^Tj 

4e 


dxi 


y/^Ti 


(53) 


Thus, 


h  < 


8k 


V^Ti 


+2Y,  [  PnAlM^i)  -  >  ^}[Wi{^i)  -  Mil  +  iMi  -  OQ\]fi{xi)dxi.  (54) 

i=i  JR 


Moreover, 


h  =  2  EE/  PnAlUxi)  -  wfe)l  >  -  Vl(^i)\ 

i=Xj=\J^i3  ^ 

X  fi{xi)fj{xj)dxidxj 

+2  E  E  / ,  PnAm^l)  - 

i=i  7=1  ^ 


X  fi{xi)fj{xj)dxidxj 

k  ft  p 

<  fi{xi)fj{xj)dxidxj 

7=1  j=l 

k  k  .  ^ 

+2  ^  /  ,  -Pn,.{|<Mi(2;i)  -  <Mi(a^i)l  > 

7=1  ^ 

X  fi{xi)fj{xj)dxidxj. 


(55) 


From  (28),  when  s  is  large  enough,  \(fi{xi)—Xi\  <  e  and  \(pj(xj)  —Xj\  <  e.  Therefore, 
when  s  is  sufficiently  large. 


{{xi,Xj)  :  \^pi{xi)  -  (pj{xj)\  <  e}  C  {(2:7,2:^)  :  \xi  -  Xj\  <  3e}.  (56) 


Thus,  similar  to  (53), 
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'f  XiA 


6e 


-\/2^inin(Ti,  Tj) 


(57) 


We  observe  that 


+  -  lJ.j\]fi{xi)fj{xj)dxidxy 


(58) 


Prom  (54)  and  (58),  it  suffices  to  analyze  the  limiting  behaviors  of 


lR^n,3{.\‘Pi{Xi)  ^i{Xi)\  >  2^ fiiX'i)dXi, 

Ir  ^Ti,s{\0iiXi)  ~  'Pi{Xi)\  >  ~  f^i\fi{Xi)dXi.  (59) 


We  first  analyze  jRPn,3{\‘Pi{xi)  -  ^pi{xi)\  >  ^}fi{xi)dxi.  Denote 


yi  =  {xi :  \(pi{xi)  -  9i\  < 

Zi  =  {xi  :  \xi  -  9i\  <  f}.  (60) 

By  Lemma  4.1,  we  know  that  when  s  is  large  enough,  \(pi{xi)  —  Xj|  <  |.  Therefore,  for 
sufficiently  large  s,  we  have 


R-yiCR-Zi 


(61) 


and 


< 


J^Pn,s{\^i{Xi)  (pi(^Xi)\  >  ^}fi{Xi)dXi 

[if  Pn,3{.\^i{.Xi)  —  ^i{Xi)\  >  fi{Xi\9i,  (7^)hi[9i\lJ,i,  T^)dXi 
JR  \JR-Zi  ^  ^ 

Lily  ^  ^}fi{xi\9i,af)hi{9i\iJ.i,Tf)dxi^ 


d9i 


d9i 
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<  fi{xi\Oh(^i)hi{Oi\[Xi,T^)dx^^dei 

+  -  ^i\  ^  2^fi{^iK<^i)hi{0M,ri)dx^^ddi 

<  f^{xi\^i,(rf)dx^^hi{^i\|Xi,Tl)dei 

^  -  Oil  2^f/(5  +  1)  >  0}fi{xi\9i,a^)dxi^ 

xhi{6i\ni,rf)d9i 

A  -  2«5-?/(s  +  1)  <  0}fi{xi\9i,  a^)dx,^hi{9i\iXi,T^)d9i 

p  (2a+l)e^ 

<2/6  128(e+l)2-*^/li(^i|^i,r/)c?^i 
^  2 


+ 


+ 


+o(n 


(2a+l)e^ 

<  2e  128(e+l)^o'*^ 


+ 


+ 


'2 

i/f  e. 


^  )  hi{9M,  ri)d9i 


+o{n  ^) 
<  0(5-1) 


+ 

+ 

+ 


Xi 

r,2 


'^i  °  '^i 

[  (f  Pn,s{%  <  i^f/(2rf)}/i(xi|^i,cr?)(ia;i]/ii(0i|//i,T2)c?0i 

J  R  \  J  R  / 


< 


xhi{9i\iJ,i,Tf)d9i 

+  ^  <  J^V{2af)}fi{xi\9i,af)dx^^hi{9i\ni,Tf)d9i 

+o(n-i) 

0(5-1)+/  f/  /i(a;i|0i,cr?)rfXi)/ii(0i|/ii,r?)d0i  +  o(n-i) 

+  /^  (/^-Pn,s{|^i  -/^i|  >  64^2^'  '^DdOi 

}fi{xi\9i,  aDdx^  hi{9i\iii,  Tf)d9i 


64cr? 

~  -  64of' 


+o(n  1) 

(23+l)./,^e^ 

<  0(5-1)  +  +  /  (,^,,.2.  hi{9i\ixi,rf)d9i 

+o{n~^)  +  o(n-i) 


64(rr 
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(62) 


(,3  +  lfvh^ 

<  0{s-^)  +  e  +  o{n-^) 

=  0{s~^)  +  o{n~^). 

Similarly,  we  can  obtain 

r  €  11 

/  PfiAlM^i)  -  >  o}lv?i(^i)  -  lii\fi{xi)dxi  =  o(-)  +  O(-).  (63) 

Jr  /  ns 

Combining  (49),  (54),  (58),  (59),  (62)  and  (63),  we  finally  obtain  the  rate  of  conver¬ 
gence  of  the  proposed  selection  procedure. 

Theorem  1.  The  selection  procedure  d(x)  defined  in  (26)  is  asymptotically  optimal 
with  convergence  rate  of  order  o(i)  -1-  0(^^).  That  is, 

B(R(S"'’Kx))  -  R(d‘(m  =  »(b  +  0(5^).  (64) 

Jh  O 

5  Simulations 

We  carried  out  a  simulation  study  to  investigate  the  preformance  of  the  selection  pro¬ 
cedure  d("'’*)(x).  The  expected  risk  £'(i?(d(”’®)(x))  —  R{d^{x)))  is  used  as  measure  of  the 
performance  of  the  selection  rule. 

We  consider  the  following  case  in  which  =  3,  that  is,  we  have  3  logistic  populations 
Hi,  112  and  Hs  and  we  would  like  to  use  the  proposed  selection  procedure  to  select  the 
best  population  compared  with  a  control. 

The  simulation  scheme  is  described  as  follows: 

(1)  For  each  n,  s  and  for  each  i  =  1, 2, 3,  generate  independent  random  variables  Xu, . . . , 
XiM  as  follows: 

'  for  /  =  1, . . . ,  n, 

<  (a)  first  generate  from  normal  distribution  with  density  A/’()ai,  rf)  (65) 
(b)  then  generate  Xiji  from  logistic  distribution  L{9ii,  cr^) 

(2)  Based  on  the  past  observations  Xiji,  and  the  present  observations  X  =  [Xi, . Xk), 
we  construct  the  empirical  Bayes  selection  procedure  d^"'’®^(x)  and  compute  the  condi¬ 
tional  difference 
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D  =  R{{d^^'^\x)  -  R(d^{x))). 


(66) 


(3)  Repeat  steps  (1)  and  (2)  400  times.  The  average  of  the  conditional  differences  on 
the  400  repetitions  which  is  denoted  by  D{n,  s),  is  used  as  an  estimator  of  the  differences 
ER{{d^^’^\x)-R{d^{x))). 

Tables  (1)  gives  the  simulation  results  on  the  performance  of  the  proposed  empirical 
Bayes  selection  procedures.  We  choose  6q  —  0.5,  =  0.4,  fj.2  —  0.5,  and  fxz  =  0.6, 

n  =  r2  =  rs  =  1. 

From  these  results,  we  see  that  D(n,  s)  decreases  to  zero  very  rapidly.  It  supports 
Theorem  1  that  the  convergence  rate  is  o(i)  +  0(^^). 


Table  1 

Performance  of  the  selection  rule 


n 

D{n,  s  =  1) 

D{n,  s  =  10) 

D{n,  s  =  50) 

5 

0.05132320 

0.01647000 

0.00560000 

10 

0.03145200 

0.00653760 

0.00218600 

15 

0.00636600 

0.00367570 

0.00079450 

20 

0.00389500 

0.00293676 

0.00010670 

30 

0.00278474 

0.00089434 

0.00008610 

40 

0.00283848 

0.00008932 

0.00004989 

50 

0.00019361 

0.00023743 

0.00003889 

60 

0.00056436 

0.00010391 

0.00004021 

70 

0.00023664 

0.00009736 

0.00002519 

80 

0.00035232 

0.06272372 

0.00001805 

90 

0.00636233 

0.00211873 

0.00001781 

100 

0.00036277 

0.00012751 

0.00001664 

125 

0.00326283 

0.00032525 

0.00001033 

150 

0.03272747 

0.00003257 

0.00000819 
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Graph  for  Table  1  (when  s  =  1) 


Silumation  (1) 
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Graph  for  Table  1  (when  s  =  50) 


Silumation  (3) 


n 
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